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ABSTRACT 

We have analyzed the long period, double-lined eclipsing binary system OGLE SMC113.3 4007 
(SC10 137844) in the SMC. The binary lies in the north-eastern part of the galaxy and consists of two 
evolved, well detached, non-active G8 giants. The orbit is eccentric with e = 0.311 and the orbital 
period is 371.6 days. Using extensive high-resolution spectroscopic and multi-color photometric data 
we have determined a true distance modulus of the system of m-M=18.83 ± 0.02 (statistical) ± 0.05 
(systematic) mag using a surface brightness - color relation for giant stars. This method is very 
insensitive to metallicity and reddening corrections and depends only very little on stellar atmosphere 
model assumptions. Additionally, we derived very accurate, at the level of l%-2%, physical parameters 
of both giant stars, particularly their masses and radii, making our results important for comparison 
with stellar evolution models. Our analysis underlines the high potential of late-type, double-lined 
detached binary systems for accurate distance determinations to nearby galaxies. 
Subject headings: binaries: eclipsing — galaxies: distances and redshifts — galaxies: individual (SMC) 
— stars: late-type 



1. INTRODUCTION 

The determination of the cosmic distance scale is one 
of the fundamental tasks in astronomy. Perhaps it is the 
most fundamental one because usually physical quanti- 
ties used in astrophysics are scaled by the distance. The 
building of the cosmic distance ladder has proven to be a 
difficult, time consuming multi-step task. The first step 
consists of using geometrical methods to measure the dis- 
tances of nearby stars and clusters. There are only a few 
well defined, purely geometrical methods we can use: 1) 
the annual trigonometric parallax, 2) the proper motions 
convergence point of a star cluster, 3) the orbital paral- 
lax of a spectroscopic visual binary star, 4) the kinematic 
parallax of a pulsar. These primary methods are cur- 
rently limited to only a few hundred parsecs around the 
sun. In some favorable circumstances this limit can be 
expanded to a few kiloparsecs: examples are the distance 
determination to the Galactic Center via the orbital par- 
allax of stars using adaptive optics e.g., iGillessen et al.l 
(2009) and the annual trigonometric parallax determina- 
tion of Milky Way mascrs using long-baselin e radio in- 
terferometry e.g.jMoellenbrock et aTT(| 2009): Rcid et all 
lj2009f) : lNagavama et al.l(|2011[ ). Toreach out beyond the 

* Based on observations obtained with the ESO NTT and 3.6m 
telescopes at La Silla, and with the Magellan Clay telescope at 
Las Campanas Observatory 



Milky Way we have to use secondary methods calibrated 
with the help of the primary ones, or use other direct but 
" near-geometrical" methods based on some astrophysical 
models. In the first case we are restricted quite severely 
by the limited volume of space where we can safely cal- 
ibrate secondary distance indicators. In the latter case 
we are affected by the model assumptions inherent to the 
different methods. Additionally, these model assump- 
tions are often based on an accepted distance scale and 
may produce a circular reasoning - see as an example the 
discussion o f the Pleiades' distance modulus problem by 
van Leeuw cn (200i|). 

A fundamental step in the construction of the dis- 
tance ladder is the determination of accurate distances 
to the Magellanic Clouds, particularly to the Large Mag- 
ellanic Cloud. This galaxy with its relatively simple ge- 
ometrical structure and smaller than Small Magellanic 
Cloud depth in the line of sight is a rich astrophysi- 
cal laboratory containing sizable samples of stellar dis- 
tance indicators of different kinds and still remains an 
ideal a nchor po i nt for the distance scale (see recent re- 
view of [Walker (2011)). One of the fundamental tasks 
of the Araucaria project is to optimize existing meth- 
ods by checking on their s ystematic errors. We have 
been using red clump s tars (Pi etrzynski fc Gieren 2002; 
Pietrzynski et al.ll2003fl . Ce pheid variables (jGieren et al.l 
2005t iStorm et al.l l2011allbf ) and RR Lyrae stars 



Graczyk, Pietrzyriski et al. 



(|Szewczvk et al.l I2008D to measure the LMC distance, 
with all these techniques taken to the near-infrared to 
minimize systematic errors due to dust extinction on the 
results. These different methods have yielded true LMC 
distance moduli in the range 18.45 - 18.58 mag, with to- 
tal errors in the range 4-6%. While this work has already 
yielded important progress in the long-standing effort of 
determining a truly accurate distance to the LMC, it 
is desirable and necessary for precision cosmology, and 
in particular for the determination of an accurate (bet- 
ter than 5%) value of the Hubble constant, to do better 
and determine the LMC distance to 1%. A very encour- 
aging step in this direction was recently made in our 
project with the analysis of a first late- type eclipsing bi- 
nary in the LMC, OGLE-051019.64-685812.3. This de- 
tached, double-lined binary system consisting of two G4 
giant stars has allowed a near-gcomctrical distance mea- 
surement accurate to 3% (jPietrzvhski et afl 12009). We 
are currently observing and analyzing a number of simi- 
lar systems in the LMC from which we expect to measure 
the LMC barycenter distance with total uncertainty bet- 
ter than 2.5%. We have also detected several detached, 
double-lined eclipsing binary systems with red giant com- 
ponents in the SMC which are extremely useful for a 
precision distance measurement to that galaxy. 

There are a large number of implementations of the 
eclipsing binary (EB) method in the literature. They 
differ mostly in amount of the modeling needed to cal- 
culate the total emergent flux of a binary system. The 
computed flux is scaled to the observed flux at Earth 
to get a distance. One implementation employs stellar 
atmosphere models to calculate the emergent flux by fit- 
ting them to the observed spectral energy distribution 
(SED) obtained with the spectrophotometry. Such an 
approach was undertaken by E. Guinan and collabora- 
tors in the middle of 90's (jGuinan et al.1 (19961 ) with the 
goal to obtain a precise distance to the LMC with a total 
uncertainty of about 2%-3%. They studied EBs with late 
O-type and early B-type components bright enough to 
obtain sufficient quality data (jGuinan et al. 1998a). Af- 
ter, their pioneering work on the distance determination 
to the early type LMC binary HV 2274 , with a claimed 
accuracy of 3.5% (|Guinan et al.lfl9~98bl ). it had become 
clear that the method has indeed a large potential. How- 
ever, the extensive use of stellar atmosphere models in 
the fitting procedure makes such an approach very model 
dependent. Because absolute fluxes, especially in the UV 
region, predicted by atmosphere models for very early 
type stars can be calibrated with an accuracy of at most 
5% (principally caused by the lack of a sufficient number 
of nearby calibrator stars with reliable geometrical par- 
allaxes) , the interpretation of the resu lts obtained from 
early type systems is uncertain, e.g. iFitzpatrick et all 
(2003 , see t h eir se ction 7) in the case of the LMC and 
iNorth et al.l (J20101 see their section 5.6) in the case of 
the SMC. 

A different possibility is to employ an empirical sur- 
face brightness - color relation and to use angular diam- 
eter scaling to obtain a distance to an eclipsing binary. 
Such relations, for different colors, are well established 
for dwarfs and giants with spectral types later than F0 
from accurate measur ements of stellar angular diame- 
ters by interferometrv (|Groenewe gen 2004; di Be nedettol 
120051: [Kervella fe Fouquel 120081; iCasagrande et al.1 12010). 



This approach was used by us to determine the 
distance to OGLE-05 1019.64-685812.3 in the LMC 
(|Pietrzvhski et al.l 12009). Although the 1% precision of 
distance determinatio n to extrag a lactic eclipsing binary 
systems envisioned by iPaczvhskl (|1997f ) is, in principle, 
possible to obtain, we are presently still hampered by the 
systematic uncertainty on the surface brightness-color 
calibration - see the section 4.1 of this paper. 

In this paper we present the distance determina- 
tion to the long period, well detached eclipsing binary 
OGLE SMC113.3 4007 (a 20 oo = 01 h 06 m lO?33; 5 20 oo = 
-72°06'25'/2) in the SMC. As this system is ideal for 
distance determination the distance, we show the full 
potential of the method and discuss all possible sources 
of error which can lower the robustness of our final re- 
sult. The star was identif ied as an eclipsing b inary in 
the catalogue published by Udal ski et al.l (J1998I ) and ob- 
tained identifier SC10 137844 in the OGLE-II database. 
Initially, because the eccentricity of the system is close 
to one third and the secondary minimum appears almost 
exactly at orbital phase 0.67, the system was incorrectly 
assu med to have a circu lar orbit and a period of ~ 248 
days. iWvithe fc Wilson! pOOll ) selected SC10 137844 as 
a prime distance indicator to the SMC and shortly af- 
terwards we started obtaining high-resolution spectro- 
scopic observations of the star. The system was iden- 
tifi ed as an eclip s ing bi nary in the MACHO database 
by iFaccioli et all {2007) as MACHO 206.17005.6. In 
the OGLE-III photometric database the star is listed as 
SMC113.3 4007. 

2. OBSERVATIONS AND DATA REDUCTION 

Johnson-Cousins optical photometry of the system was 
obtained with the Warsaw 1.3 m telescope at Las Cam- 
panas Observatory in the course of t he second, third and 
fourth phases of the OGLE project (jUdalski et al.lll997t 
lUdalskill2003l) and also with the ANDICAM camera at- 
tached to the CTIO 1.3 m telescope. In total we se- 
cured 1178 I-band epochs and 168 V-band epochs. Be- 
cause of the long orbital period consecutive epochs were 
taken usually on different nights. The time span of the 
I-band observations is 5119 days (from JD 2450627 to JD 
2455746). The raw d ata were reduced with the im age- 
subtraction technique (|Wozniakll2000l: lUdalskill2003t) and 
instrumental magnitudes were calibrated onto the stan- 
dard system using Landolt standards. We did not use 
MACHO data for this star because OGLE I-band pho- 
tometry has better eclipse coverage and better precision 
(smaller noise and much smaller number of outliers) . 

Near-infrared photometry was collected with the ESO 
NTT telescope on La Silla equipped with the SOFI cam- 
era. The setup of the instrument, and the reduction 
and calibration of the data onto the UKIRT system 
were essentially identical to those descr ibed in our pre- 
vious paper by iPietrzvhski et "all (|2009l ). The transfor- 
mation of our photometry onto th e Johnson s ystem was 
done using the equatio ns given bv lCarpenterl ()200lD and 
IBessell fc Brettl (p88h . 

High resolution echelle spectra were collected with the 
Clay 6.5 m telescope at Las Campanas Observatory, 
equipped with the MIKE spectrograph, and with the 
3.6 m telescope on ESO La Silla, equipped with the 
HARPS spectrograph. We employed an 5 x 0.7 arc sec slit 
with MIKE giving a spectral resolution of about 40000. 
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Table 1 

Radial Velocity Measurements for SC10 137844 



4884 4886 

Wavelength (A) 



4890 



Figure 1. A MIKE spectrum of OGLE SC10 137844 taken at the 
third quadrature. Some absorption lines from both components are 
identified. 

In the case of HARPS we used the EGGS mode yield- 
ing a resolution of about 80000. In our analysis we used 
8 HARPS spectra, 16 MIKE blue region and 14 MIKE 
red region spectra (25 epochs; 13 times blue and red 
MIKE spectra were taken simultaneously). The typical 
S/N at -5000 A was about 20 and 10 for MIKE and 
HARPS spectra, respectively. Figure Q] shows a small 
part of a MIKE blue region spectrum taken close to the 
third quadrature (JD 2454816.6). 

Before further analysis all the spectra were shifted to 
the Solar System Barycenter. In order to derive ra- 
dial velocities of the system's component s we employed 
the Broadening Function (BF) formalism (Rucinski 1992, 
119991 ) and the two-dimensional cross-c orrelation method 
TODCOR bv lMazehfc Zuckerl (U99l . We used numer- 
ous metallic lines in the wavelength regions 4125-4230, 
4245-4320, 4350-4840, 4880-5000, 5350-5850, 5920-6250, 
6300-6390, 6600-6800 A. As templa tes we used syntheti c 
spectra from a library computed bv lCoelho et all (|2005| ). 
The templates were chosen to closely match the atmo- 
spheric properties of the stars in a grid of T a s and logg. 
We iterated the radial velocity determination a few times 
each time updating the template according to the latest 
solution. Usually the difference between velocities de- 
rived by TODCOR and BF was smaller than 100 m s" 1 
and we adopted the mean value of both as final mea- 
surement. We detected a linear trend in velocities de- 
rived from MIKE spectra taken between HJD 2454000 
and 2454800 equal to about —0.6 m s _1 per day. After 
removing the trend we accounted also, during the men- 
tioned period of time, for a zero point shift between the 
MIKE and HARPS spectrographs of 0.344 km s" 1 . The 
individual RV measurements are listed in Table [TJ 

3. MODEL OF THE SYSTEM 

To derive absolute fundamental parameters for the sys- 
tem we used two codes for the analysis of the eclipsing 
binary sta rs: the Wilson-Devinney (WD) program, ver- 
sion 2007 dWilson fc Devinnevl H97ll: 1 Wilson! fl979L fl990l: 
Ivan Hamme fc Wilson! I2007T ) equipped with the auto 
mated Differential Correction (DC) optimiz ing subrou- 
tine; and 
2004aB) 



the JKTEBOP program dSouthworth et al 
based on the EBOP code 
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1981!) . The WD code allows to simultaneously solve for 



* HARPS spectra 
+ MIKE red spectra 



radial velocity and multi-band light curves recommended 
as the best way to obtain a consistent model of a binary 
system, e.g.. iWilsonl (J1979D . while use of JKTEBOP al- 
lows for a reliable determination of the errors by Monte 
Carlo simulations. The WD code is based on Roche-lobe 
geometry and employs quite a sophisticated treatment 
of stellar surface physics, while in the JKTEBOP code 
stellar surfaces are approximated by biaxial ellipsoids, 
and simpler description of stellar physics is used. As our 
system is well detached it is meaningful to compare the 
results from both programs. 

3.1. Parameter Choice 

The choice of adjustable parameters has crucial im- 
pact on the final solution and parameters of the model. 
We tried to estimate and fix as many parameters as 
possible. The preliminary orbital period was calculated 
with the string-length method giving P = 371.71 days. 
The moment of the primary minimum was estimated to 
be T = 2450776.684 from the I-band light curve and 
fixed during the later analysis. In this paper we will 
refer to the primary as the star which is eclipsed in 
the deeper, primary minimum. The average tempera- 
ture of the primary was estimated as follows. We de- 
termined the reddening in the direction of our target to 
be E(B-V)= 0.06 ± 0.03 usin g the OGLE-II redden- 
ing maps (lUdalski et alJll999t) and OGLE-III redden- 
ing maps (Haschke et al.ll2011h . Average out-of-eclipse 
magnitudes were calculated from all observations taken 
outside minima: F(15?807), /(14™722), J(14?020), 
if(13™317). These magnitudes were der eddened using 
the in terstellar extinction law given by ISchlegel et al.l 
(1998) and assuming Ry = 3.1. From combined prelim- 
inary solutions obtained from RaVeSpAn software writ- 
ten by B. Pilecki and from JKTEBOP code we obtained 
the surface gravities: loggi = 1.7 and log #2 = 1.6, the 
eccentricity e = 0.31 and T^ — T\ ~ 10 K. The metallic- 
ity was assumed to be [Fe/H]= —0.5 appropriate for a 
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Table 2 

Effective Temperatures Derived from the Dereddened Colors 



Table 3 

Photometric Parameters 



T.//(JO 



Calibration 



Color Index 



4732 di Benedetto (1998) V-K 

4763 Alonso et a~ (7999) V-K 

4799 Houdashelt et al"72000') V-K 

4717 Ramirez k, MelendezT2005') V-I,V-K 

4830 Masana et al. (2006) V-K 

4847 Gonzalez H. fc Bon ifacio (2009) V-K 

4899 Casagrande et al. (201CT) V-I,V-K 

4812 Worthev fc Lee (20lTT J-K,V-I,V-K 

4803 mean 



you ng (age ~ 200 My r) population of stars in the SMC: 
e.g. iDias et all (J2010T ) from integral spectroscopy of clus- 
ters for a cluster NGC 458 (it lies very close to our star) 
derived [Fe/H]= -0.4 and age about 100-200 Myr (sim- 
ilar to age of our star). To derive an average effective 
temperature of the binary we employed a number of em- 
pirical calibrations published during last thirteen years 
(Table [5]) . The resulting average effective temperature 
was set as a temperature of the primary Ti = 4800 ± 100 
K. It corresponds to a sp ectral type G8 giant star accord- 
ing to the calibration by I Alonso et al.l (|1999f ). The lower 
metallicity e.g., [Fe/H]= —0.7 has only minor effect on 
the temperature (it is lower by 10 K). 

Using BF we estimated the velocity broadening to be 
8.4 km s" 1 and 10.8 km s _1 for the primary and the sec- 
ondary, respectively, with errors ±1.0 km s^ 1 . These ve- 
locities are lower than the appropriate semi-synchronous 
velocities expected if the components are in tidal lock- 
ing during periastron passage, 12.5 km s _1 and 13.2 km 
s _1 , respectively. If we account for the macroturbulence 
velocity field typical for a G ty pe giant atm osphere de- 
scribed by Crt = 5.4 km s" 1 (|G7avll200l Table B.2), 
the primary velocity broadening is consistent with syn- 
chronous rotation: vs'mi = 6.4 km s _1 . However, the 
secondary seems to rotate super-synchronously with the 
rotational velocity vsmi being by a factor of 1.4 larger 
than the synchronous one. Accordingly we set the star 
rotation parameters F\ — 1.0 and F-i = 1.4. 

The albedo parameter was set to 0.5 and the gravity 
brightening to 0.32, both values appropriate for a cool, 
convective atmosphere. The limb darkening coefficients 
were calculated internal ly by the WD c ode according 
to the logarithmic law of IKlinglesmith feS obicski (1970) 
during ea ch iteration o f DC using tabulated data com- 
puted by Ivan Hammd (|1993l) . With the help of BF we 
tried to detect any signature of a tertiary component in 
the spectra at different orbital phases but we failed, and 
thus set the third light parameter to I3 = 0. 

As free parameters of the WD model we chose the semi- 
major axis a, the orbital eccentricity e, the argument of 
periastron u, the phase shift of the primary spectroscopic 
conjunction cj>, the systemic radial velocity 7, the orbital 
inclination i , the secondary star average surface temper- 
ature T2, the modified surface potential of both compo- 
nents ill, ^2, the mass ratio q = M2/M1, the observed 
orbital period P b s , and the relative monochromatic lu- 
minosity of the primary star in the two bands Lly, LI/. 
It is worth noticing that within the WD code there is 
no possibility to directly adjust the radial velocity semi- 
amplitudes K\ and K%, instead the semimajor axis and 



Parameter 



WD a 



JKTEBOP b 



Eccentricity e 

Argument of Per. u) (deg) 

Phase Shift <j> 

Orbital Inclination i (deg) 

Sec. Temperature T% (K) 

Fractional Radius n 

Fractional Radius T2 

(rl + rl) 

k = r2/rl 

Observed Period P i, 3 (d) 

(L2/L1) V 

(L2/L1)! 

(L2/Ll)j 

(L2/Ll) K 

Third Light I 3 



0.3112 ±0.0012 

30.2 ±0.4 

0.07557 ± 0.00009 

88.19 ±0.02 

4813 ±3 

0.1097 ±0.0007 

0.1160 ± 0.0007 

0.2257 

1.057 

371.768 ± 0.002 

1.135 ±0.011 

1.131 ±0.012 

1.125 c 

1.120 c 

0.0 



0.3064 ± 0.0010 

28.5 ±0.5 

not adjusted 

88.18 ±0.02 

not adjusted 

0.1145 ±0.0016 

0.1117 ±0.0016 

0.2262 ±0.0004 

0.975 ±0.028 

371.768 ± 0.002 

not adjusted 

0.976 ±0.052 



not adjusted 



a Simultaneous solution of V-band and I-band light curves, and ra- 
dial velocity curves of both componets. Uncertainties quoted are the 
standard errors from DC subroutine. 



Only I-band light curve solution. 
Monte Carlo simulations. 
c Extrapolated from the model. 



Uncertainties quoted are from 



the mass ratio are adjusted simultaneously. 

3.2. Fitting Procedure 

We fitted simultaneously two light curves, in the I-band 
and V-band, and two radial velocity curves using the DC 
subroutine of the WD code. The detached configuration 
(Mode 2) was chosen during all the analyses and a sim- 
ple reflection treatment (MREF=1, NREF=1) was em- 
ployed. A stellar atmosphere formulation was selected 
to both stars (IFAT=1). A level dependent weighting 
was applied (NOISE=l) and curve dependent weightings 
(SIGMA) were calculated after each iteration. The ini- 
tial input solution was found using the JKTEBOP code. 
Convergence was defined to have been achieved if the pa- 
rameter corrections given by DC were smaller than their 
standard errors on three consecutive iterations. The grid 
size was initially set to N=40 for both stars but other 
grids were checked. It was found that very fine grids 
(N> 50) produced some numerical instability and no con- 
vergence, in the sense defined above, could be obtained. 
We did not break the adjustable parameter set up into 
subsets, but instead we adjusted all free parameters at 
each iteration. 

At the end of the fitting procedure we additionally ad- 
justed the third light I3 to find its impact on the solu- 
tion: the third light corrections were invariably negative 
suggesting an unphysical solution, and we therefore set 
I3 = in our final solution. 

The I-band light curve solution obtained with the WD 
code is presented in Figures [2] and [3] This solution was 
compared with the solution resulting from the applica- 
tion of the JKTEBOP code. Using JKTEBOP we fitted 
only the I-band light curve, and values of the mass ratio 
and the limb darkening coefficients were set according 
to the final WD solution. The errors of the parameters 
were calculated from ten thousand Monte Carlo simula- 
tions (Mode 8 of JKTEBOP code). It turned out that the 
parameters obtained from both solutions are consistent 
with each other within 2cr errors with one exception: the 
individual fractional radii of the components. The WD 
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Table 4 

Orbital Parameters 



Parameter 



Value 



Periastron Passage (HJD) 2451114.80 

Orbital Period P (d) 371.585 ± 0.002 a 

Semimajor Axis a (R Q ) 417.49 ± 0.79 

Systemic Velocity 7 (km s _1 ) 147.68 ± 0.05 

Prim. Velocity Semiamplitude K\ (km s" 1 ) 29.93 ± 0.10 

Sec. Velocity Semiamplitude K 2 (km s _1 ) 29.46 ± 0.11 

Mass Ratio q 1.016 ± 0.005 

1 Corrected for the radial movement of the centre of the system mass 
n respect to Solar System Baryccntcr - see Scction l3.3l 
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Figure 2. The I-band light curve solution to OGLE SC10 137844 
from the WD code. 
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Figure 3. Zoom into the I-band light curve solution to OGLE 
SC10 137844. Left panel presents the primary minimum, the right 
panel shows the secondary minimum. 



solution prefers the secondary star being the larger and 
more luminous of the two stars, while the JKTEBOP so- 
lution prefers the opposite. It underlines the old problem 
of determining univocally stellar radii for eclipsing bina- 
ries having partial eclipses. To distinguish between the 
solutions we investigated integrated CCF and BF power 
profiles of both stars to derive their spectroscopic light 
ratio. The calculated secondary to primary light ratio 
from CCF and BF turned out to be about 1.1. Thus, 
for consistency of the model, we decided to adopt the 
radii delivered by the WD code but with uncertainties 




0.4 0.6 
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Figure 4. The V-band light curve solution to OGLE SC10 137844 
from the WD code. 
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Orbital phase 

Figure 5. The radial velocity curve solution to OGLE SC10 
137844 from the WD code. Filled black circles denote the primary's 
velocities. The dashed line marks the systemic velocity. The devi- 
ation of the secondary radial velocities from the Keplerian motion 
near the orbital phase 0.69 indicates the Rossiter-McLaughlin ef- 
fect. 



as returned from the JKTEBOP code. The resulting 
ephemeris for the primary minimum is: 

T (HJD)=2450776.68(±0.01) + £x371.7684(±0.0017) 

The V-band solution is shown in Figure HI and the 
radial velocity curve solution in Figure [SJ The photo- 
metric and orbital parameters of the final solution are 
given in Table [3] and Table |H respectively. We note 
that the systemic velocity of our system is in excellent 
agreement with the mea n radial velocity of t he red giant 
branch stars derived by lHarris fc Zaritskvl (|2006[ ) from 
their IMACS field 3. As that field encompasses the eclips- 
ing binary, the systemic velocity independently confirms 
that our target belongs to the SMC. 

Two of our spectra were taken during secondary eclipse 
egress and have allowed an approximate modeling of the 
Rossiter-McLaughlin effect. The observed moderate pro- 
file (see Figure [SJ was fitted by adjusting the secondary 
star rotation parameter and we obtained F2 = 1.3 ± 0.1. 
This value seems to confirm that the secondary rotates 
super-synchronously and it is in accordance with the ini- 
tial value F2 = 1.4 derived from the BF analysis. 
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Table 5 

Physical Properties of the OGLE SC10 137844 
System 



Property a,b 


The Primary 


The Secondary 


Spectral Type 


G8 II-III 


G8 II-III 


V (mag) 


16.630 


16.493 


V — I (mag) 


1.087 


1.083 


V-K (mag) 


2.498 


2.483 


J — K (mag) 


0.705 


0.701 


Radius (Rq) 


45.8 ±0.7 


48.4 ±0.7 


Mass (M ) 


3.504 ± 0.028 


3.561 ±0.025 


logff (cgs) 


1.660 ±0.017 


1.620 ±0.016 


T eff (K) 


4800 ± 100 


4813 ± 100 


vs'mi (km s 1 ) 


6.4 ±1.0 


8.7 ±1.2 


Luminosity (Lq) 


1000 ± 80 


1130 ±80 


M bo l (mag) 


-2.75 


-2.88 


My (mag) 


-2.40 


-2.54 


Fe/H 


-0.5 dex 


(assumed) 


E(B-V) 


0.06 ±0.03 



a Absolute dimensions were calculated assuming: G — 
6.673 ■ 10~ s cm 3 g _1 s" 2 , R G = 6.9551 • 10 10 cm, M = 
1.9888 ■ 10 33 g, M holtQ = +4.75. 

The magnitudes and colors are observed values. 



3.3. Absolute Dimensions 

The observed orbital period of the system P b s and the 
true orbital period P are linked through a relation: 



Po, 



bs 



P(l+1) 



(1) 



where c is the velocity of light. The corrected period P 
is listed in Table 0] and is used to calculate semimajor 
axis of the system. Although the correction is small it is 
much larger than our precision of period determination. 
Once the semimajor axis a of the system is known it is 
possible to obtain the absolute dimensions of the stars. 
Table [5] gives astrophysical data about the two compo- 
nents. The observed individual magnitudes and colors 
were calculated from the out-of-eclipse mean magnitudes 
and the light ratios given in Table [3] The position of 
the disentangled components on the color-magnitude di- 
agram (CMD) is shown in Figure [6] The physical radii 
of the stars result from the relation: R = r ■ a, where r 
is the fractional radius listed in Table [3l The masses are 
derived from the equations: 



Mi[M ] = 1.34068 • 10" 

M 2 [M @ ] = Mi -q 



1 [Rq 



1 + q P 2 [d] 



(2) 
(3) 



Derived absolute visual luminosities of both componets 
are consisten t with luminosity class II-III for G-type star 
(Eggen 1994) and their masses lie in the mass range ex- 
pected for bright giant stars which have masses from 
about 3 M to 9 M Q . 

4. DISTANCE TO THE SYSTEM 

To derive the distance of our binary system we em- 
ployed a surf ace brightness - color relation. We used the 
calibration of Idi Benedetto! (|2005l ) between F-band sur- 
face brightness and V — K color relation (SBR) obtained 
from a number of carefully selected mixed sample of giant 
and dwarf stars having their angular diameters measured 
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Figure 6. The observed color-magnitude diagram of the OGLE 
field SMC113.3 containing the eclipsing binary SC 10 137844. Su- 
perimposed (full circles) are the positions of its two giant compo- 
nents. Both components of our system are clump giants located 
at the red loop of the helium burning. At this quiet and relatively 
long lasting evolutionary stage stars do not show significant pho- 
tometric and spectroscopic variations, which make them very good 
targets for precise determination of their stellar parameters. 



Table 6 

Distance determinations to OGLE SC10 137844 from different 
surface brightness-color relations 



d [kpc] 



Luminosity Class Color 



Reference 



57.5 ± 2.0 

58.0 ±2.8 

61.6 ±4.2 
60.8 ±2.0 


Dwarfs+Subgiants 

Giants 

Dwarfs 
Dwarfs+Subgiants 


V- 

V- 
V 

V- 


-K 
-K 
-I 
-K 


58.4±1.8 


Giants 


V- 


-K 



Kervella et al. (.2004) 

Groenewegcn (2004) 

Kervella &: Fouq uc (2008) 

Casagrande et al. (2010) 

di Benedetto (2005) 



precisely by long baseline interferometry. The angular di- 
ameter of a star predicted by the surface brightness-color 
relation takes the form: 



:1 Q0.2-(S-m ) 



(4) 



where S is the surface brightness in a given band and 
nio is the unreddened magnitude of a given star in this 
band. The distance in parsecs follows then directly from 
angular diameter scaling and is given by a simple linear 
equation: 



d[pc] =9.2984 



R[Rq] 
<b [masl 



(5) 



The distance to the OGLE SC10 137844 system was cal- 
culated as the average distance to both components, and 
its value is 58.4 ± 1.8 kpc corresponding to a true dis- 
tance modulus of 18.83 ± 0.06 mag. The distance moduli 
of the primary and the secondary differ by just 0.001 
mag, serving as an independent check on the model con- 
sistency. We compare our result with distance estimates 
obtained from other recent surface brightness - color cal- 
ibrations, results being summarized in Table [6) Within 
the errors all distances are consis tent with our adopt ed 
value based on the calibration bv ldi Benedett o (2005). 

4.1. Error Budget 
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The systematic uncertainty on the surface brightness 
relation itself provides the largest cont ribution to the to- 
tal err or of 2.0%, (compare with e.g., iPietrzvhski et all 
(2009)). The metallicity dependence of the SBR is very 
weak, for o ur system we expect a correction of about 
0.006 mag {Thompson et al.ll200H: fdi Benedetto! [l998h . 
and we decided to add an additional error of 0.3% in the 
error budget. Systematic zero point uncertainties of ab- 
solute photometry are 0.7% for near infrare d SOFI bands 
and 0.5% for the optical OGLE V band (jUdalski et al.l 
2000). The estimated total systematic error is 2.4%. 

The absolute radii are known with an accuracy of 1.5%. 
It would mean that the contribution to the distance error 
from the radius uncertainty was 1.5%. However, because 
of the similar diameters and temperatures of the compo- 
nents this contribution is much smaller and, to the first 
order of magnitude, equal to the uncertainty of the sum 
of the radii. Following Equation [5] the average distance 
to the system can be written as: 



Table 7 

Comparison of recent distance determinations to the SMC. 



(d>=| 



Rk, 



R 2 



Ri 



±R2 + RiAf\ {6) 
h 4>2 <i>\ J 



where index 1,2 denotes, the primary and the secondary 
star, respectively, A(f> = (f)2 — <t>i ancl k = 9.2984. The 
second term in Equation [6] constitutes, in our case, only 
about 3% of the first term. Regarding the error contri- 
bution from the radius uncertainties the first term gives 
^{r 1+ r 2 )/{R\ + #2) = 0.35/94.3 i.e. about 0.4% and the 
second term gives 0.05%. Combining it with the semima- 
jor axis uncertainty of 0.2% results in a total uncertainty 
of 0.5%. 

The accuracy of V — K color determination is about 
0.015 mag adding 0.7% to the statistical error. The dis- 
entangling of individual V-magnitudes and near infrared 
color indexes contributes only slightly to the total error 
due to the similar temperatures of the stars. We recal- 
culated the distance modulus of the system by enlarging 
the light ratios in the V and K bands with the error 
quoted in Table [3] and an error of 0.014, respectively. 
The resulting distance moduli differ by just 0.0005 mag. 
Thus we omitted this uncertainty in the error budget. 

In the fitting procedure we set atmosphere model for- 
mulation in the WD code. We have checked how the 
distance estimate will change if we assume a blackbody 
approximation. We re-fitted the data setting IFAT=0 
(black body) to both stars. The resulting radii and light 
ratios lead to a value of the distance modulus which dif- 
fers by less than 0.001 mag from our adopted value. 

The interstellar absorption, because of the reddening 
vector being almost parallel to the SBR, contributes rel- 
atively little to the total error: 0.7%. A different inter- 
stellar extinction law with Ry = 2.7 contributes insignif- 
icantly on the level of 0.3%. Adding the contributions 
of these errors on the distance determination to OGLE 
SC10 137844 quadratically we obtain a total statistical 
error of 1.1%. 

We conclude that the uncertainty of the extinction and 
V — K color determination are the main contributors to 
the statistical error in the distance. The calculation of 
the total uncertainty is a bit ambiguous because it de- 
pends on the systematic uncertainty probability distribu- 
tion which we don't know a priori. We can give however 
the limits for the total error. The lower error limit of 



(m - M) 


Stat. Err. 


Sys. Err. 


Method 


Ref a 


18.91 


0.04 


~0.1 


Eclipsing Binaries 


1 


18.96 


0.10 


? 


RR Lyr 


2 


18.93 


0.02 


? 


Bump Cepheids 


3 


18.97 


0.03 


0.12 


RR Lyr 


1 


18.85 


0.07 


~0.1 


Type-II Cepheids 


5 


19.11 


0.03 


1 


Eclipsing Binaries 


6 


18.92 


0.14 


1 


Classical Cepheids 


7 


18.83 


0.02 


0.05 


Eclipsing Binary 


8 
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2.5% comes from the assumption of a gaussian distribu- 
tion of the systematic uncertainty while the upper limit 
of 3.5% results from the assumption of a pure bias. The 
mean value of both was given in Table |5] corresponding 
to the total uncertainty of ~3.0%. 

5. DISCUSSION 

We obtain a true distance modulus for our binary 
system of 18.83 ± 0.02 (statistical) ± 0.05 (systematic 
from the SBR) mag. Our analysis demonstrates that for 
such late-type systems consisting of red giant stars we 
can indeed measure the distances very accurately, with 
prospects of improving the achievable accuracy once an 
improved surface brightness-color relation becomes avail- 
able which is the current limiting factor on the precision 
of the method. 

In Table 7 we have compiled other recent determina- 
tions of the distance to the SMC. Our result from SC10 
137844 is in good agreement with the average SMC dis- 
tances derived from classical and Type-II Cepheid, and 
RR Lyrae star sample s. The outlier is t he distance of 
19.11 mag reported by INorth et al.l (|2010f) from an anal- 
ysis of 33 eclipsing binaries in the southwest part of the 
galaxy. However, since the SMC is a line-of-sight elon- 
gated structure, especially towards t he northeast, with 
a typi cal depth of ~ 0.25 mag (e.g. lHarris fc Zaritskvl 
(2006)), all the determinations in Table 7 might be com- 
patible within their quoted errors. Also it is known 
that there is a distance gradient over the SMC with its 
northeastern pa rt being closer to us than its southwes t- 
crn part (e.g. ISubramanian fc Subram aniaml ([2011)). 
Indeed, the distanc e to the galaxy centre given by 
IHilditch et al.l (120051) li es between our determination and 
the INorth et al.l (|2010l ) estimate which seems to be con- 
sistent with the pr esent geometrical model of the galaxy. 
However, because INorth et al.l {2010) did not give any 
estimate of the systematic uncertainty of their result it 
is not easy to compare the SMC distances obtained from 
eclipsing binaries. 

If we assume a distan ce gradient over the SMC pro- 
jected disc of 2-3 kpc (Hildit ch et al.l [2005) our star 
should be around ~1 kpc (0.035 mag) closer to us than 
the optical center of the galaxy. Thus our distance would 
correspond to a distance modulus of 18.87 mag for the 
SMC center. Assuming a modulus offset between SMC 
and LMC equal to 0.44 (jCioni et al.ll2000D it would imply 
a distance modulus of the LMC of 18.43 mag. This value 
is very close to the true LMC distance modulus derived 
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recently by iStorm et al.l (|2011bl ) with the near-infrared 
Baade-Wesselink method applied to classical Cepheids. 

We are currently working on a number of similar late 
type eclipsing systems in the Small Magellanic Cloud, 
and expect to determine the mean distance to the SMC 
with a precision better than 4%. 
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